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We propose a new type of a three-loop induced neutrino mass model with dark matter 
candidates which are required for the neutrino mass generation. The smallness of neutrino 
masses can be naturally explained without introducing super heavy particles, namely, much 
heavier than a TeV scale and quite small couplings as compared to the gauge couplings. We 
find that as a bonus, the anomaly of the muon anomalous magnetic moment can simultane¬ 
ously be explained by loop effects of new particles. In our model, there are doubly charged 
scalar bosons and leptons from isospin doublet fields which give characteristic collider signa¬ 
tures. In particular, the doubly charged scalar bosons can decay into the same-sign dilepton 
with its chirality of both right-handed or left- and right-handed. This can be a smoking 
gun signature to identify our model and be useful to distinguish other models with doubly 
charged scalar bosons at collider experiments. 

PACS numbers: 


I. INTRODUCTION 

One of the most important questions in the standard model (SM) for elementary particle physics 
is why neutrinos have non-zero but quite small masses as compared to those of other fermions. This 
question can be interpreted to ask why the coefficient of the dimension five Weinberg operator [lj] 
(c/A)LiLi<f >< h (c and A being dimensionless and dimensionful parameters, respectively) is quite 
small, which gives Majorana neutrino masses after replacing the Higgs field <f> by the vacuum 
expectation value (VEV). From the various measurements of neutrino oscillations, typically (c/A) 
must be the order of 10~ 14 GeV -1 . 

The seesaw mechanism fsj] gives us the simple explanation why we have such a small coefficient, 
namely, a right-handed neutrino is quite heavy, e.g., 0(1O 14 ) GeV with a 0( 1) Yukawa coupling. 
However, to test this mechanism, we need a super high energy collider, and it seems to be unrealistic. 

As an alternative scenario, there are radiative neutrino mass models originated from the Zee 
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TABLE I: The particle contents and charge assignments. The superscript i denotes the flavor index (i = 1-3). 


model |3|. In these models, the smallness of the coupling c can be naturally explained by a loop 
suppression factor l/(167r 2 ) 7V in a IV-loop model. As a result, we do not need to introduce super 
heavy particles such as the above mentioned right-handed neutrino, and thus we can test these 
models at collider experiments. In addition, in radiative neutrino mass models, some of new 
particles which appear in a loop diagram for the neutrino mass generation can be a dark matter 
(DM) candidate as in the model by Ma |d ]. Empirically, three loop models can provide appropriate 
size of the coefficient c with TeV scale new particles and order of 0.1-1 couplings. The model by 
Krauss, Nasri and Trodden Q] is the first proposed three loop model with DM. In addition to those 
models, it has been proposed other three loop models in Refs, jfl-l]]. 

In this Letter, we build a new type of a three loop neutrino mass model with DM candidates. 
To have three loop diagrams, we introduce hypercharge Y = 3/2 isospin doublet scalar and lepton 
fields which include doubly and singly charged components. As a bonus, the anomaly in the muon 
g — 2 can also be simultaneously explained by loop effects of the Y = 3/2 particles. Characteristic 
signatures at collider experiments are expected from the decay of the doubly charged scalar bosons, 
where they can mainly decay into the same-sign dilepton with the both right-handed and left- and 
right-handed chirality. This signature can be used to distinguish the other models with doubly 
charged scalar bosons. 


II. THE MODEL 

We discuss a three loop neutrino mass model in the electroweak SU(2)l x U(1)y gauge theory 
with an additional unbroken discrete Z 2 symmetry. The particle contents are shown in Table Q] 
We extend the lepton sector from the SM, i.e., we introduce vector-like lepton fields Ly 2 with 
Y = —3/2 and right-handed neutrinos (i = 1-3). The scalar sector is composed of isospin 

doublet fields with Y = 3/2 ($ 3 / 2 ) and with Y = 1/2 (<h), and two singlets S + with Y = 1 and S° 
with Y = 0. For simplicity, we take S° to be a real held. 
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The Zi symmetry plays three important roles in our model. First, it forbids the L R e&*N R term 
(e being the 2 by 2 anti-symmetric matrix) which provides the tree level Dirac neutrino mass term. 
Second, it also forbids the T 3 / 2 el ^* e R term which gives the tree level mixing between the exotic 
charged leptons and the SM charged-leptons. Finally, the Z 2 symmetry guarantees the stability of 
the lightest neutral Z 2 odd particle, i.e., N R or S' 0 , so that it can be a candidate of DM. 

In this setup, the lepton sector Lagrangian is given by 

-Aep = \m' n N%N' r + Mi /2 {L\ /2 ) L (U Z/2 ) R + h.c. 

+ VSM L L® e R + Vr e R N R S+ + h ' c - 

+ fi r L^ /2 ) R S + + f$ + h.c., (1) 

where the first and second terms in the first line correspond to the Majorana masses for N l R and 
the Dirac masses for Ly 2 , respectively. In the second line, y SM represents the usual Yukawa 
matrix which provides the masses of the SM charged-leptons. All the other terms give new Yukawa 
interactions which are necessary to obtain left-handed Majorana neutrino masses excepted for the 
y R term. 

The most general Higgs potential is given by 

Hffiggs = + P 3 / 2 1^3/21 2 + P±\S 2 + 2^o(5'°) 2 

+ Ail^l 1 + A | <1*3/2 1 4 + AalS+l 1 + A 4 ( 5 0 ) 4 

+pii^i^i 2 + p-2\®\ 2 \ s+ \ 2 + P m 2 (sy 

+ P4|^ > 3/2| 2 |‘S' + h + P51^3/2 PCS' 0 )" + /?6|‘S' + | 2 (S' 0 )" 

+ ^‘M 2 + [^l / 2 ^S + S° + &l /2 e$( 5 -) 2 + h.c.], (2) 

where the scalar doublet fields <!> and ^3/2 can be parameterized as 


$ = 

G + 

, ^*3/2 — 

1- 

+ CO 

-e- 

1_ 


- j=( v + h + iG°) 


1 

-e- 

to 

1_ 


with v (~ 246 GeV) being the VEV, and G± and G° being the Nambu-Goldstone bosons which 
are absorbed into the longitudinal components of and Z bosons, respectively. In our model, 
h corresponds to the discovered Higgs boson with the mass of about 125 GeV. We note that the 
singly-charged scalar bosons (S^, (j>y 2 ) do not mix with each other because of the Z 2 parity. Among 
the parameters in the potential, n and £ terms play a crucial role for the neutrino mass generation 
and the decay of ( h^ 2 , respectively, as we shall see below. 
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FIG. 1: Three-loop neutrino mass diagrams. The red lines denote the Z 2 odd particles. 


The squared masses m 2 of physical scalar bosons tp(= < h^ 2 , S' 1 * 1 , S°, h ) are then simply 

given without introducing mixing angles as follows 

m 2 h = 2 An 2 , 

u 2 v 2 

m |±± = M3/2 + m |± = P3/2 + ir(Pi + (J )) 

3/2 Z 3/2 ^ 

v 2 

m l± = *4 + yP2, m 2 s o = l4 + v 2 P3- (4) 

We here note that in our model, only the Higgs doublet field <L gets the non-zero VEV, so that the 
masses of the SM weak gauge bosons as well as charged fermions are given exactly the same form 
as those in the SM. 


III. DARK MATTER 


The lightest neutral Z 2 odd particle can be a DM candidate; i.e., N l R or S' 0 . Let us first suppose 
Nft to be the lightest particle among the Z 2 odd particles. In this case, only the t- and u-channel 
processes i + i~ being the SM charged lepton) contributes to the DM annihilation 

via the mediation at the tree level. This gives the p-wave amplitude, and the cross section is 
calculated by 


<7U re l 


\Vr\ A 4 + x 2 2 
487 tM 2 n1 (1 + x) 4 reb 


( 5 ) 


where x = m 2 s± /M‘^ 1 and n re i is the relative velocity of DM. We note that there is no s-wave 
contribution due to the chiral suppression. To reproduce the observed relic density, the cross 
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section given in Eq. (JSJ) should be inside the following region 


av re i = (1.78-1.97) x 1(T 9 GeV 


T -2 


( 6 ) 


at the 2a level. This can be satisfied by taking y T \ to be 0(1 ) with the mass of DM of 0(0.1-1) 
TeV. 

Next, we give a brief comment on the bosonic DM ( S °) case. This corresponds to the so-called 
Higgs portal DM scenario, where DM is annihilated via the s-channel Higgs boson h mediation. 
In Ref. [loj], it has been found that the DM mass of about 500 GeV can satisfies the relic density 
without conflicting the constraint from the direct searches such as the LUX experiment [11]. In 
addition to this solution, when the DM mass is around the half of the Higgs boson h mass, i.e., 
~ 63 GeV, both the relic density cross section and the direct search experiments can be satisfied 
using the h pole effect. 


IV. NEUTRINO MASS 


In our model setup shown in Table [T] the tree level contribution to the left-handed neutrino 
masses, i.e., via the type-I seesaw mechanism is forbidden by the Z 2 symmetry. In addition, 
there are no one- and two-loop diagrams, where a systematic classification for one- and two-loop 
diagrams for Majorana neutrino masses has been done in Refs. [12, 13]. As a result, the leading 
order contribution to the neutrino masses is given at the three loop level as depicted in Fig. [TJ 
where the topology of these diagrams is the same as neutrino mass diagrams in Ref. Qj. 

There are four types of contributions to the neutrino masses which are denoted by the diagrams 
(I), (II), (III) and (IV). The contributions of (I) and (II) pick up two Dirac masses of L 3 / 2 , while 
those of (III) and (IV) do not. We note that the diagram (I) gives a different contribution from 
the diagram (II), because of the difference of the loop momentum flow. Similarly, diagrams (III) 
and (IV) provide a different results with each other. 

The mass matrix for the neutrinos is then separately calculated by the four parts: 


(■— 
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k 2 v 2 


(16tt 2 ) 3 2M 2 


X - 
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(M l X c + (M 1 ^ + (M l } l )$ c + (M™) 


{ II\a6c 
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v )ij 


I IV \abc 
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(7) 


where M max denotes the largest value of the masses among , N^, S± and S°. (A4]f v )fj c 

correspond to the contributions from the diagram I-IV, respectively. Each contribution is calculated 
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= ^ac/i°(/i)“ 6 M^(4) bc /f Fjy, 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 


where Epyi denote the loop functions for each diagram. Although the loop functions are determined 
by fixing all the mass parameters of the particles running in the loop diagram, their typical values, 
i.e., 0 ( 0 . 1 ) are not so sensitive to these mass parameters. 

The neutrino mass matrix is parametrized by 


= (CIpmns 0pMNs)*i> 


( 12 ) 


where m]) ias = diag(mi, m 2 , m 3 ) is the diagonalized neutrino mass matrix with three mass eigen¬ 
values mi, m 2 and n 73. The Pontecorvo-Maki-Nakagawa-Sakata matrix [l4, 15] C/pmns is the 3x3 
mixing matrix to diagonalize the neutrino mass matrix. A global fit on the recent neutrino os¬ 
cillation data Q] provides their two squared mass differences and three mixing angles of I/pmnSj 
depending on the mass ordering. The total neutrino mass Ylh is constrained from the experi¬ 


ment by Planck Collaboration 


17], which suggests m k < 0.23 eV at 95 % CL. We expect that 


observed three mixing angles can easily be reproduced by the Yukawa couplings Jl, f r and /lr 
as seen in Eqs. (| 8 |) - (fTT|) . In particular, the case with diagonal /V- and fn couplings is favored to 
avoid lepton flavor violating processes, and even in such a case, the mixing can be accommodated 
by /lr- 

We here estimate typical values of the parameters that appear in the above neutrino mass 
formulae. Assuming that each of Al[; TV is the order of 0.1 eV as required by the magnitude of the 
observed neutrino masses, we obtain the following relation for the contribution, e.g., from 


1 v 2 M\ /2 M n 

(T^F * F ' ' = ° (tu) eV ’ 


(13) 


where the flavor indices are omitted here. When we take M max ~ M 3 / 2 — IWjv, and F I ~ 0.1, 
Eq. (1131) can be rewritten as 


(«/ l/lrY 


= 0( 10 " 5 ). 


Mr, 


(14) 
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FIG. 2: New contributions to the muon g — 2. In the left diagram, the photon can be emitted from either 
L ~/- or S~. 


In the case of M max = 0(1) TeV, the factor (h/l/lr) 2 should be the order of O(10 -4 ) which can 
be realized by taking each of k. fj J and Jlr couplings to be in the range of 0.1 to 1. Therefore, 
tiny neutrino masses can be naturally explained by introducing neither super heavy particles nor 
quite small coupling constants as compared to the gauge coupling constants. 

V. MUON ANOMALOUS MAGNETIC MOMENT 


In our model, it is possible to explain the anomaly in the muon anomalous magnetic moment 
(muon g — 2) as a bonus. It has been known that there is a discrepancy between the SM prediction 


and the measured value at Brookhaven National Laboratory. The observed value is given by 


m 


a e * p = 11659208.0(6.3) x 10 


-10 


and the difference between the above value and the SM prediction (A p — a)) M ) was calcu¬ 
lated as 


Aa M = (29.0 ±9.0) x 10 _1 ° 
A= (33.5 ±8.2) x 10 -10 


19 ] 

20 ] 


and 


(15) 

(16) 


The above results given in Eqs. (fl5ll and (fl6l) correspond to 3.2<r and 4.1cr deviations, respectively. 
There are new contributions to Aas shown in Fig. O These are calculated as 

Aa M = ^| l/ ^ |2 [ G ( M 3/2’ m s ± )+ 2G ( m 5±,^)] -\y^G{M^m s± ) 1, (17) 


where 


G(mi,m 2 ) =- 


2 m\ ± 3mfm2 — 6± — 6mfm2 In 


6 (m\ — m^) 4 


(18) 
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M 3/ / [GeV] M 3/ / [GeV] 


FIG. 3: Parameter region which can explain the discrepancy of the muon g — 2 on the My 2 -m s ± plane. 
The black, blue and red shaded regions give the prediction of Awithin the lcr, 2a and 3er error given in 
Eq. (USD. The value of /^ is fixed to be 1.5, 2.0, 2.5 and 3.0 in the upper-left, upper-right, lower-left and 
lower-right panels, respectively. In all the plots, we take y^ = 0. 


In this loop function, we neglect the muon mass dependence. In Eq. (1171) . we take both fi and yn 
couplings to be a diagonal form in order to avoid constraints from lepton flavor violating processes 
such as y —> ey 1 . Thus, only L and contribute to the muon g — 2. 

In Fig. [3] we show the parameter region which can explain the anomaly of the muon g — 2 on 
the Myr,-m s ± plane. We here use the discrepancy of Agiven in Eq. (fT6l) . In these plots, we take 
yff = 0 to remove the contribution from which always gives the distractive effect to the 
loop contribution as seen in Eq. (H71) . In this case, Ais given as a function of m s ±, M.^ 2 and 
f 1 ^ 1 . The value of /^ is fixed to be 1.5, 2.0, 2.5 and 3.0 in the upper-left, upper-right, lower-left 
and lower-right panels, respectively. The parameter region indicated by black, blue and red shaded 


1 There are £ —> 31 type lepton flavor violating processes via box diagrams [2lj, j22|. However, these contributions 
can also be avoided by taking diagonal couplings of ya and /l . 
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regions gives the prediction of Awithin the la, 2 a and 3a error. We find that the 2a allowed 
region is obtained by taking rn s ± + < 400, 600, 750 and 950 GeV for /^ = 1.5, 2.0, 2.5 and 

3.0, respectively. 


VI. COLLIDER PHENOMENOLOGY 


In this section, we discuss collider signature of the new particles in our model. The most 
important signature to identify our model could come from the doubly charged scalar bosons 
and leptons L^. At the LHC, these particles are produced in pair via the Drell-Yan process, i.e., 
pp — > ( I l 'jy '2 ^ 3/2 an ^ PP ^ ^ 3/2 ^. 3/2 ’ an d a ^ so v ^ a assoc i a ted processes with singly-charged particles, 
i.e., pp -)• and pp -)• Lf±L* /T 

The decay pattern of these doubly charged particles strongly depends on the mass spectrum of 
the new particles which can be roughly determined from the discussion in the previous sections and 
experimental constraints. From the discussion in Sec. V, we found that we need to take masses of 
^ 3 / 2 ^ anc ^ to a f ew hundred GeV to explain the anomaly in the muon g — 2. In addition, 
if we consider the case of the scalar DM, i.e., S°, m s o should be taken to be about 63 GeV as 
we explained in Sec. III. Although in this case the masses of N\ i are not constrained so much, it 
should be at most 0(1) TeV to reproduce the order of neutrino masses as we discussed in Sec. IV. 
For the mass of the most important bound comes from the same-sign dilepton searches at 

the LHC. The current bound using the LHC Run-I data gives the lower limit on the mass of 
to be about 550 GeV when mainly decay into the same-sign dielectron [zsj]. Weaker bounds 
are obtained if the flavor of the same-sign dilepton from the decay of is different from the 
electron. We will discuss later how the same-sign dilepton decay of is realized in our model. 
Regarding to the doubly charged leptons L^, we cannot apply the similar bound for the mass of 
*^ 3/2 that of Ly 2 , because their decay products must include the DM in the final state due to 
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those Z 2 odd property. Thus, the typical final state of the decay of L^p is the same-sign dilepton 
plus missing energy, and the mass bound must be weaker than that of Tijyjp The study for the 
collider phenomenology of doubly charged leptons have been studied in Ref. 24]. 

Under this configuration, let us consider the following mass spectrum: 

M n (~ 0(1) TeV) > m $3/2 > M 3/2 > m s ± > m s 0 (~ 63GeV), 


(19) 


where = m$±±(= m $ ± ) 2 , and we take the flavor universal masses M 3 / 2 = M^p 2 ’ 3 and 


3/2 


3/2 


Mn = M'' n 1,2,3 for simplicity. In this mass spectrum, the following decay modes can be considered: 


<F 


±± 


3/2 

T ±± 
-^ 3/2 


S ± S ± (via £), <h ^ 2 —> (loop induced), 


( 20 ) 

S ± e$ (via f L ), ( 21 ) 

with the subsequent decay 3 of S± &±* 2 S° —> via the k, term in the potential. In (1201) . the 

loop induced decay into the same-sign dilepton is driven via the diagram shown in Fig. |4j This 
process can also be understood by the following dimension five operators with a typical mass scale 
of our model M: 


\ l l $ 3/2 d^e R , 


' e R e R^3/2^, 


( 22 ) 


M 0/ " M 

where the first (second) operator induces the left (right) decay process of Fig. |4l The important 
point here is that can decay into the same-sign dilepton with its chirality of both right-handed 
or left- and right-handed. This feature cannot be seen in other doubly charged scalar bosons from 
the other models as follows. First, doubly charged scalar bosons from an isospin singlet with Y = 2, 
which are introduced in the two loop neutrino mass model the so-called Zee-Babu model { 26 ], can 
decay into the same-sign dilepton with only right-handed. Next, doubly charged scalar bosons 
from an isospin triplet with Y = 1, which is introduced in the Higgs triplet model (HTM) j^, can 
decay into the same-sign dilepton with only left-handed. Therefore, by measuring the chirality of 
the same-sign dilepton event, our model can be distinguished from the other models with doubly 
charged scalar bosons. In Ref. [28(, it has been clarified that the same-sign dilepton signals with a 


different chirality can be distinguished by using distributions of decay products of tau leptons. 

Finally, we comment on the signal and background events in our model. If the same-sign 
dilepton decay of is dominant which can be realized by a smaller value of the £ coupling in 


2 If there is a non-zero mass difference between m and m ± , the cascade decay, i.e., <I>’T —> W is 

*f/2 *f/2 3/2 3/2 

possible. In Ref. [25(, the collider phenomenology for such a cascade decaying doubly charged scalar bosons has 
been discussed. 

3 The —> Iv decay is given via similar diagrams as shown in Fig. [4] 
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the potential, the expected collider signature is pp —>• Z* j 7 * —> ^ 3/2 ^ 3/2 ^ • Dedicated 

studies of the signal and background analysis for the four lepton events has been done in Refs. 


29- 


311 ]. The SM background processes ZZ, WZ, tt, Zbb, Ztt and Wtt have been taken into account 


in Ref. 


311 ] in the detector level analysis. Among them, the largest background cross section comes 


from the ZZ and Ztt processes, but they can be significantly reduced by taking the invariant mass 
cut of the i + t~ system M^+£~, i.e., the event with \M(+(- — m z \ < 10 GeV is excluded [3l|]. As 
a result of the background separation from the signal, the 5a discovery reach for doubly charged 


Higgs bosons in the HTM has been presented in [3l|, where the required integrated luminosity is 
6 fir -1 (60 ftu 1 ) when the mass of doubly charged Higgs bosons are taken to be 500 (700) GeV 
at the 14 TeV LHC. In our model, the production cross section of is smaller than that in 
the HTM because of the difference of the doubly charged Higgs boson coupling to the Z boson. 
For example, if we take the mass of doubly charged Higgs bosons to be 500 (700) GeV, we obtain 


the pair production cross section of 1 (0.15) fb and 2 (0.3) fb in our model and in the HTM 


m, 


311 ] can be applied in our 


respectively. However, the same signal and background analysis given in 
case, so that we expect to obtain the 5a discovery reach for by requiring the roughly 2 times 
larger integrated luminosity than the case in the HTM with the same mass of the doubly charged 
Higgs bosons (due to about the 50% reduction of the cross section as compared to the HTM) at 
the 14 TeV LHC. 


VII. CONCLUSIONS 


We have constructed a three-loop induced neutrino mass model with the additional unbroken 
discrete Z 2 symmetry. The lightest neutral Z 2 odd particle, i.e., the right-handed neutrino N }j or 
the real singlet scalar boson S° can be a DM candidate. We have shown that the neutrino masses 
can be naturally explained by the order 1 TeV new particles with the order 0.1-1 couplings due to 
the three-loop suppression. At the same time, the anomaly of the muon g — 2 can also be explained 
by the loop effect of the doubly charged leptons . Regarding to the collider phenomenology, 
we have discussed the signature from doubly charged scalar bosons 4 >% 2 and leptons In 

particular, we have emphasized that the loop induced decay of into the same-sign dilepton 
an d can important signature to distinguish our model from the other models 

with doubly charged scalar bosons. 
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